Time Structure of Ultra-High Energy Cosmic Ray Sources and Consequences for 

Multi-messenger Signatures 
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The latest results on the sky distribution of ultra-high energy cosmic ray sources have consequences 
for their nature and time structure. If the sources accelerate predominantly nuclei of atomic number 
A and charge Z and emit continuously, their luminosity in cosmic rays above ~ 6 x 10 19 eV can be 
no more than a fraction of ~ 5 x 10~ 4 Z~ 2 of their total power output. Such sources could produce 
a diffuse neutrino flux that gives rise to several events per year in neutrino telescopes of km 3 size. 
Continuously emitting sources should be easily visible in photons below ~ 100 GeV, but not in TeV 
7— rays which are likely absorbed within the source. For episodic sources that are beamed by a 
Lorentz factor T, the bursts or flares have to last at least ~ 0.1 T -4 A~ 4 yr. A considerable fraction 
of the flare luminosity could go into highest energy cosmic rays, in which case the rate of flares 
per source has to be less than ~ 5 x 10~ 3 F 4 A 4 Z 2 yr -1 . Episodic sources should have detectable 
variability both at GLAST and TeV energies, but neutrino fluxes may be hard to detect. 
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I. INTRODUCTION 

The sources of ultra- high energy cosmic rays (UHECR) 
above the "GZK threshold" Q around ~ 6 x 10 19 eV are 
still unknown. Recently, an important step forward has 
been made by the Pierre Auger Observatory which has 
revealed a correlation of the arrival directions of UHECR 
above ~6x 10 19 eV with the nearby cosmological large 
scale structure as mapped out by the distribution of ac- 
tive galactic nuclei (AGNs) [2]. At least if UHECR de- 
flection in large scale cosmic magnetic fields is moderate, 
this requires a certain minimal density of sources within 
the "GZK horizon" of about 75 Mpc. At the same time 
the observed spectrum normalizes the required injection 
power per volume. Together, these two numbers imply an 
upper limit on the time averaged UHECR injection power 
per source. Comparing this with the minimal power that 
needs to be dissipated in order to produce UHECR up to 
~ 10 20 eV, this allows to constrain the time structure of 
the UHECR emission, in particular continuously emit- 
ting versus episodic sources. This can also have some 
implications on how the dissipated total power may be 
distributed between the cosmic ray, photon, and neutrino 
channels. The latter can be important for a multimes- 
senger study of UHECR sources. In the present paper we 
attempt to work out these constraints in a largely model- 
independent way, including their dependence on the type 
of nuclei that are predominantly accelerated. 

The remainder of this paper is structured as follows. In 
Sect. II we develop general requirements on the individual 
sources. In Sect. Ill and IV we consider continuously 
emitting and episodic sources, respectively. In Sect. V 



we discuss Centaurus A as a potential UHECR source 
and we conclude in Sect. VI. We will use the units in 
which c = 1 throughout. 



II. REQUIREMENTS ON INDIVIDUAL 
SOURCES 



Accelerating particles of charge eZ to an energy E max 
requires an induction £ > E max /(eZ). With Zq — 100 
the vacuum impedance, this requires dissipation of a min- 
imal power of [!, 0| 
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L m \n — TT- — 10 Z 



V 10 20 eV 



ergs 
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We stress that this minimal power can be less by factors 
of order unity in specific geometrical circumstances, such 
as magnetic fields connecting accretion disks with jets 
111 AGNs @. However, given other, larger uncertainties 
such as the chemical composition of UHECRs, we can 
ignore such details in the present work. 

The "Poynting" luminosity Eq. |T]) can also be ob- 
tained from the expression L ni i n ~ T 2 (BR) 2 where T 
is the beaming factor of the accelerating region and the 
product of the size R and magnetic field strength B of the 
acceleration region is given by the "Hillas criterium" Q 
which states that the Larmor radius tl = E max / (TeZ B) 
should be smaller than R, 
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(2) 
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In the following we denote cosmic rays above 6 x 10 19 eV 
as ultra-high energy cosmic rays (UHECR) and we take 
£" max ~ 10 20 eV as the benchmark for their typical pro- 
duction energy within the sources. Any source producing 
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UHECR up to energy E max at a given time has to have 
a total power output of at least the Poynting luminosity 
Eq. ([1]). Note that this is comparable to the Eddington 
luminosity L Edd {M) = 1.3 x 10 38 (M/M Q ) ergs -1 of a 
massive black hole of mass M in the centers of active 
galaxies. A considerable part L 7 of that power is pre- 
sumably electromagnetic and thus emitted in photons. 
We now assume that electromagnetic power is produced 
in the same area of size R in which UHECR are accel- 
erated. Denoting the characteristic photon energy by s, 
the optical depth for pion production on such photons by 
accelerated protons with an energy above the photo-pion 
threshold, £>6.8x 10 16 (e/eV) _1 eV, is given by 
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where we have used <r P7 ~ 300/Aarn around the thresh- 
old for pion production. Note that R ~ 1 pc is the 
typical size of an accretion disk around a supermas- 
sive black holes at the centers of AGNs which is de- 
termined by the "sphere of influence" ~ 2G-^M/v 2 ~ 
2(M/10 7 M Q )(w s /200kms- 1 )- 2 pc where G N is New- 
ton's constant and v s is the velocity dispersion of the 
stars in the host galaxy 0] • The optical depth for photo- 
disintegration of primary nuclei is comparable to Eq. ([3]). 
If it is significantly larger than unity, most nuclei would 
be disintegrated before leaving the source and the max- 
imal energy would have to be -E max ^ A 10 20 eV in or- 
der for UHECR to arrive at Earth with energies up to 
10 20 eV. 

On the other hand, the optical depth for hadronic in- 
teractions of accelerated protons and nuclei with the sur- 
rounding bulk matter of hadronic mass Mb u ik extending 
over a characteristic scale i?t> u ik ~ R can be written as 



^ CppMbuik 
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where we have estimated the nucleon density by n p ~ 
-^buik/-Rb u ik- Note that the mass of AGN accretion disks 
is roughly comparable to the mass of the central super- 
massive black hole Q whose typical mass is 1O 7_8 M . 
Eq. is only a rough estimate because the details will 
depend on the geometry, for example spherical versus 
disk-like accretion. Since the bolometric luminosities of 
most AGNs are < 10 47 ergs 1 , a comparison of Eqs. ([3]) 
and ((4]) suggests that hadronic interactions dominate over 
photo-hadronic interactions in the cores of AGNs. 

Pionic and photo-hadronic processes will produce sec- 
ondary 7— rays and neutrinos. The optical depth for 
photons of energy above the pair production threshold, 
E > m 2 Je ~ 0.26(e/eV)- 1 TeV, can be estimated as 
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where or — 0.6 barn is the Thomson cross section. 

We now deduce some requirements on the size R of 
the accelerating region. We will also take into account 
a possible beaming factor T such that B and R and 
other length scales are measured in the comoving frame, 
whereas luminosities and the energy E max refer to the ob- 
server frame. The synchrotron loss length for a nucleus 
of atomic number A and charge Z in a magnetic field of 
strength B is 



^synch 



0.43 r A z~ 
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V 10 20 eV 



B 

g' pc - 



(6) 



and the Larmor radius can be written as 



r L ~0 



IT^Z- 1 ( -4 
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V10 20 eV/ \G 



pc. (7) 



Since any energy loss time must be longer than the accel- 
eration time which itself is larger than the Larmor radius, 
one has the condition ^ S y nc h ^ which gives an upper 
limit on the magnetic field strength, 



B < A.3T 2 A A Z 
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V.10 20 eV 
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Together with Eq. {7J this results in a lower limit on the 
Larmor radius, 



> 2.3 x IQ-^^A 



E n 



V 10 20 eV 



pc. (9) 



In case of AGN sources, for example, this is certainly 
consistent with a size R ~ 1 pc for the typical size of an 
accretion disk. Acceleration could thus occur in small 
parts of such accretion disks. 



III. CONTINUOUSLY EMITTING SOURCES 

Assuming at most moderate deflection in intergalactic 
space, the number of arrival directions observed by the 
Pierre Auger Observatory and other experiments implies 
a lower limit on the source density 0, Q , 



> 3 x 10 _5 Mpc" 



(10) 



The UHECR flux observed by the Pierre Auger observa- 
tory is flpj ] 



<1Nck 
dE 



(E ~ 6xl0 19 eV) ~ 6xlO~ 40 cm" 2 sr" 1 s" 1 eV _1 , 

(11) 



which corresponds to a power per volume of Q 

Quhe ~ 1.3 X 10 37 erg Mpc~ 3 s _1 . (12) 
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Comparing Eqs. (|T0|) and ([12]) implies for the time- 
averaged UHECR luminosity per source 



J XJHE 



< 4 x 10 41 ergs" 1 



(13) 



This is much smaller than the instantaneous total lumi- 
nosity required by Eq. {!]). 

If UHECR sources emit continuously, Eqs. ([I]) and (fl3|) 
imply that these sources must emit at least ~ 2000 Z~ 2 
times more energy in channels other than UHECR. This 
is consistent with the fact that at redshift zero an average 
AGN in an active state has a bolometric luminosity of 
~ 5 x 10 44 erg s _1 , comparable to Eq. |T]), and the volume 
emissivity is ~ 3 x 10 40 erg M pc~ 3 s _1 , a factor of a few 
thousand larger than Eq. |12|) [Tllll2l|. The average AGN 
bolometric luminosity and volume emissivity corresponds 
to a density of "typical" AGNs of ~ 6 x 10- 5 Mpc~ 3 , 
consistent with Eq. (fTTT)) . 

As a result, if sources emit continuously and the total 
power is distributed roughly equally between hadronic 
cosmic rays and electromagnetic power, the cosmic ray 
injection spectrum could extend down to < 10 17 eV with 
a rather steep spectrum oc E~ a , a < 2.7. Using Eq. (TT2")) , 
this can be written as 



dn C R 
dE 



(E)~(a-2) 



E 



QuHE 

(1 20 eV )2 yiQSOeV 



(14) 



Furthermore, Eq. (0| suggests that the optical depth for 
hadronic interactions can be of order unity and thus a 
considerable part of that cosmic ray flux could be trans- 
formed to neutrinos with energies ~ 10 17 eV. Following 
Ref. [l3l |. for proton primaries, Z = 1, we can write for 
the production rate per volume of neutrinos 



dhv it?\ 2 / dn CR 
-dE (E ^3^^E- iE/Xv) > 



(15) 



where x v ~ 0.05 is the average neutrino energy in units of 
the parent cosmic ray energy and / = e T — 1 is the ratio 
of number of cosmic rays interacting within the source to 
cosmic rays leaving the source. If the cosmic ray injection 
spectrum oc E~ a extends down to -Emm without break, 
/ is limited by 



V 10 20 eV 



< 2 x 10 3 Z- 2 



(16) 



where L to t is the total luminosity and L m ; n is given by 
Eq. (fTJ). This condition just results from comparing the 
total emissivity with the output in UHECR and neutri- 
nos and would be saturated if the total output would 
be dominated by neutrinos in case of "hidden sources". 
Since neutrinos do not interact during propagation and 
ignoring redshift evolution, we can estimate the all-flavor 
diffuse neutrino flux as 



1 dh v 
4irH ~dE 



(E) 



(17) 



where Hq = 100/ikms 1 Mpc 1 is the Hubble constant 
with h ~ 0.72. Putting together Eqs. (QUI, US]) and 1(17). 
we obtain 



E 2 j*> s (E) ~ 190<" 1 (a-2)/ 

2— a 



„a — 1 ^ 

E 

10 20 eV 



(18) 



eV cm 2 s 1 sr 1 



If the ankle marks the transition from galactic to extra- 
galactic cosmic rays, the injection spectral index of the 
latter has to be a ~ 2.2, especially if heavier nuclei are 
accelerated 14| . The secondary neutrino spectrum would 
then extend down to at least ~ 10 17 eV and Eq. (fl~8|) im- 
plies 



E 2 jf*{E)~1.2f 



E 



10 17 eV 



-0.2 



eV cm 2 s 1 sr 1 



(19) 

In contrast, if the ankle is due to pair production of ex- 
tragalactic protons, then one needs a ~ 2.6 15]. The 
secondary neutrino spectrum would then extend down 
to at least ~ 10 16 eV and Eq. JTHJ) implies 



E 2 j« m (E) - 240/ 



E 



10 16 eV 



-0.6 



eV cm 2 s 1 sr 1 



We now compute the neutrino event rates in kilome- 
ter scale neutrino observatories for these two scenar- 
ios. Using the neutrino-nucleon cross section a U N 



1.9 x 10- 33 (£/10 16 eV) 



\0.363 . 



for 10 16 eV < E < 



10 21 eV [16|, we obtain the rate 

Ru ^ <J uN {E)2nE 3 ^{E)n N V cS 



(21) 



2.3 



E 



-0.637 



10 16 eVy Vkm 3 
E 2 jf s (E) 



100eVcm _2 sr x s 1 



yr 



where rijv — 6 x 10 23 cm -3 is the nucleon den- 
sity in water/ice and V e s the effective detection vol- 
ume. The scenario a — 2.2 with neutrino flux down 
to 10 17 eV would give ~ 2 x 10~ 2 / yr" 1 km -3 < 
20 Z~ 2 (L to t/ L min ) yr" 1 km -3 , where we have used 
Eq. (fT6|) for /. The scenario a — 2.6 with neutrino 
flux down to 10 16 eV would give ~ 5.5 f yr^ 1 kmT 3 < 
180Z- 2 (L tot /L min )yr" 1 km- 3 . 

TeV 7— rays would mostly be absorbed within the 
source due to Eq. ^ . As opposed to Ref. [l3| , we there- 
fore do not get a constraint from the non-observation 
of AGNs at TeV energies in this scenario. In contrast, 
X-rays and GeV 7— rays could leave the source. Indi- 
vidual sources should be visible by X-ray telescopes and 
by GLAST. In fact, EGRET has seen a diffuse flux [l7| 
which constrains the neutrino flux because a compara- 
ble amount of energy goes into photons and neutrinos in 
primary cosmic ray interactions, 



E 2 fi ltt (E) < 10 3 eVcm~V 



(22) 
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Eqs. (|19|) and (|20[) satisfy this limit except for sources 
deeply in the hidden regime, / ^> 1. 



IV. FLARING SOURCES 

If the sources flare on a typical time scale St in the 
observer frame, the corresponding life time of the burst 
in the comoving frame, TSt must be larger than the co- 
moving acceleration time scale which itself is larger than 
the Larmor radius, thus with Eq. ([9]) we have 



St > 0.1T- A A- 



I- 



V 10 20 eV 



yr. 



(23) 



This is consistent with the variabilities observed for 
AGNs which are observed at time scales down to ~ 
60s provided that r > 20 and/or predominantly 
heavier nuclei are accelerated. The time scale is also con- 
sistent with 7— ray bursts which can easily have Lorentz 
factors r > 20 

If we denote the rate and UHECR luminosity of typical 
flares by Rf and £uhe > respectively, we can write for the 
time averaged UHECR power, Eq. (Tie 



RfStL 



UHE 



^UHE 



< 4 x 10 41 ergs 



(24) 



This means that the fraction of time a typical intermit- 
tent source emits the typical instantaneous UHECR lu- 
minosity iuHE, also called duty factor, is given by 



V = R t St<4xlO- 4 Z 2 (I^l] (J£?- 



(25) 



Eqs. (f2"4"]l and (|25p imply that a considerable fraction 
of the minimal flaring luminosity Eq. (fT|) can go into 
UHECR, Luhe ~ imin, which could be the case if the 
UHECR acceleration spectrum is hard, a < 2. From 
Eqs. ([M]) and we then have 



Rf <5xl0- 3 T 4 A 4 z 2 



\L V heJ VlO 20 eV 



yr 



(26) 

In the limit Luhe — > ^uhe we obviously recover the limit 
of continuous sources, RfSt — > 1. 

During one flare the total non-thermal energy release 
would be 



Et > L m - m 5t > 3 x 10 51 r~ 4 A- 4 



E u 



10 20 eV 



erg, 
(27) 

If this energy release is due to accretion onto a 
central black hole with an energy extraction effi- 
ciency of ~ 10% [13], this corresponds to about 
0.02 T -4 Z~ 2 M Q . In AGN scenarios involving ac- 
cretion, this energy requirement is certainly modest. 

Individual sources would be observed with the appar- 
ent UHECR luminosity 



^UHE.obs — -E'UHE 



St 

delay 



(28) 



where idciay ^ Z 2 10 5 yr is the time delay of charged cos- 
mic rays due to deflection in cosmic magnetic fields. 

TeV 7— rays may or may not be observable from in- 
dividual sources because the duty cycle is small and 
most of the time the source would have luminosi- 
ties < 10 45 Z~ 2 ergs -1 , leading to fluxes < 2 x 
10 _8 Z -2 (d/20Mpc)" 2 ergcm- 2 s- 1 where d is the dis- 
tance to the source. In the active phases, Eq. ([5]) suggests 
that TeV 7— ray may be absorbed by pair production 
within the sources. Note that flares in the electromag- 
netic luminosity are not expected to correlate with the 
UHECR luminosities due to the large UHECR time de- 
lays. 

In the flaring limit we can have Luhe ~ -^min in which 
case both 7— ray fluxes and secondary neutrino fluxes can 
not be much larger than the UHECR flux which would 
then be a considerable fraction of the total energy bud- 
get. Since the spectrum must be rather hard in this case, 
both the diffuse and discrete neutrino fluxes are likely 
unobservably small. 



V. CENTAURUS A AND OTHER AGN 
SOURCES 

Centaurus A is the nearest AGN at a distance d ~ 
4Mpc with a central supermassive black hole of mass 
M ~ 10 8 M Q [2l|. It is not a blazar as its jet has a 
large inclination angle to the line of sight. The Pierre 
Auger Observatory measured two UHECR events from 
the direction of Centaurus A [2j. This corresponds to a 
flux [H 



dE 



(fi~6x 10 19 eV) < 10~ 40 cm 



•- 1 s- 1 eV- 1 , 
(29) 

and to an apparent UHECR luminosity of Centaurus A 
of iuHE.obs ^ 10 39 ergs _1 . If Cen A is a continuous 
UHECR source, this is consistent with Eq. (JT3j) . The 
bolometric luminosity of Cen A is L to t — 10 44 ergs _1 
which originates mostly within ~ 500 pc from the center 
and is mostly emitted around energies e ~ 1 eV [2l|, [23| . 
This is consistent with Eq. (TT|) if predominantly heavier 
nuclei are accelerated, Z > 4. At MeV energies Cen A 
has a luminosity ~ 10 42 ergs _1 [2l| . 

If Cen A is emitting UHECR continuously, and as- 
suming the UHECR injection spectrum is cx E~ a , the 
expression analog to Eq. (Ti"5"|) for a discrete source gives 
for the secondary neutrino flux 



ME) 



2/ dN CR 
3xv dE 



(E/o 



(30) 



Using Eq. (|29|) . one obtains numerically for the neutrino 
energy flux 



E 2 j v {E)~Q.2Ax a v - 1 f 



E 



6 x 10 19 eV 



,t -2 -1 

ev cm s 



(31) 
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For a ~ 2.2 this gives 



^OE) ~ 0.024/ 



E 



10 17 eV 



eVcm^s" 1 , (32) 



where from Eq. (|16|) . / <■ 100. For a ~ 2.6 it yields 



E 2 j v {E) -0.37/ 



10 16 eV 



-0.6 



eVcm^s" 1 , (33) 



where from Eq. (fT75)) . / < 3.2. The latter, more 
optimistic case would give an event rate of ~ 1.3 x 
lO^/yr^knU 3 < 4.2 x 10~ 3 yr" 1 km -3 . It is clear 
that the rate due to the diffuse flux, estimated below 
Eq. (f2"5| . is always much larger as long as CeN A is an 
" average" source. This is consistent with the conclusion 
in Ref. H. 

If Cen A is an episodic UHECR source, as may be sug- 
gested by its variability on time scales of days observed 
in X-rays and 7— rays [21[, Eqs. |28|) and ([23]) imply for 
the UHECR luminosity during a flare, 



L 



UHE 



< 



10 45 T A A 



A Z 2 



^delay \ / 



EL 



Z 2 10 5 yr / V 10 20 eV 



ergs 
(34) 

Comparing with Eq. ([T]) this suggests that the UHECR 
flare luminosity is comparable to the total output, as 
long as the flare duration is not much larger than the 
theoretical minimal variability time scale Eq. (I23[) . 

The closest blazars whose jets are close to the line of 
sight and thus may have considerably beamed emission 
are in general too far away to be responsible for UHECR. 
As an example we briefly discuss Markarian 501. This 
blazar at a distance d ~ 130 Mpc shows emission up to 
TeV energies with variability on time scales of days and 

46, 



peak luminosities of close to 10 ergs 



25] |26j. This 

is consistent with Eq. fl} even for protons, Z = 1. The 
power of such blazars is, therefore, certainly sufficient 
to provide the UHECRs. Even if they produce UHECR 
only in flares, the flare luminosity in UHECR, Luhe, 
could be a small fraction of the total flare luminosity, 
and the necessary flaring time scale Eq. (|23|) and rate 
Eq. I|28p would be consistent with observations, especially 
for significant beaming factors T typical for blazars. 



VI. CONCLUSIONS 

We have discussed some consequences of latest results 
on ultra-high energy cosmic rays for the nature and vari- 



ability of the sources as well as for the secondary 7— ray 
and neutrino fluxes produced within the sources. To this 
end we assumed predominant acceleration of nuclei of 
atomic mass A and charge Z. In the limit of continu- 
ously emitting sources, their luminosity in cosmic rays 
above ~ 6 x 10 19 eV can be no more than a fraction of 
~ 5 x 10~ 4 Z~ 2 of the total source power. If these cos- 
mic rays are produced in the accretion disks in the centers 
of AGNs, significant neutrino fluxes could be produced 
by hadronic interactions, especially in scenarios in which 
extragalactic protons dominate down to ~ 10 17 eV such 
that the ankle is due to pair production of these protons. 
The resulting cosmological diffuse neutrino flux can lead 
to detection rates up to several events per year and km 3 
of effective detection volume. This also implies consider- 
able photon fluxes at energies up to — 100 GeV, the latter 
of which should be easily visible by GLAST. In contrast, 
TeV 7— rays are likely absorbed within the source. For 
episodic sources that are beamed by a Lorentz factor T, 
individual flares have to last at least ~ 0.1 T~ 4 A~ 4 yr. 
Such flares can also be visible in photons up to the TeV 
energy range. A considerable fraction of the flare lumi- 
nosity could go into highest energy cosmic rays which 
suggests a hard injection spectrum. In this case the rate 
of flares per source has to be < 5 x 10~ 3 T 4 A 4 Z 2 yr -1 . In 
contrast to continuously emitting sources, both neutrino 
fluxes from individual sources and the resulting cosmo- 
logical diffuse flux may be hard to detect in the limit 
of flaring sources. Conversely, if high energy neutrinos 
are soon detected, this may suggest sources that produce 
ultra-high energy cosmic rays continuously. 
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